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^ ■ Abstract 

, A theory of nucleus-nucleus collisions has been developed for kinetic energies sub- 

stantially in excess of the binding energy. The very high pressure produced in the 
" . compound system as a result of the fusion of the two colliding nuclei is the reason for 

I the subsequent hydrodynamic expansion of the nuclear medium. The energy and an- 

gular distributions of the reaction products are investigated. The charge distribution 
■ is also determined in the case where the nucleon and ion components of the reaction 

products are predominant. A solution is found for the expansion into vacuum of a 
sphere in which the initially uniformly distributed material is initially at rest and at 
an ultrarelativistic temperature. 
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^ ■ 1 Introduction 
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Progress in the technology of acceleration of multiply charged ions [1] has substantially 
contributed to recent developments in this important field of research in nuclear physics. 
The complexity of the colliding systems, i.e., the accelerated ion and the target nucleus, 
gives rise to a variety of possible reaction channels specific for this category of processes. 
Let El be the kinetic energy of the incident nucleus per nucleon. For E'l ~ 1 -h 10 MeV, the 
nucleon binding energy in the initial systems, and the Coulomb barrier, which impedes the 
approach of the two particles, may still play an appreciable role. This, of course, leads to an 
increase in the fraction of reactions involving the transfer or capture of individual nucleons 
during the interaction [2]. However, these values do not, in principle, represent the limit of 
experimental possibilities, and there is a promising tendency for Ei and the atomic weight 
of the colliding systems to increase. The physical picture may then be expected to undergo a 
substantial change, and the predominant mechanism responsible for most of the interaction 
cross section turns out to be relatively simple. 

To avoid unnecessary detail with very little bearing on the essence of the situation, we 
shall confine our attention to a head-on collision between two identical nuclei and, unless 
stated to the contrary, we shall carry out our analysis in the center-of-mass system (c.m.s.). 

*Zh. Eksp. Teor. Fiz. 70, 768-784 (1976) [Sov. Phys. JETP 43, No. 3, 397-405 (1976)] 
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Suppose that the kinetic energy per nucleon in this system is Eq, the mass number of the 
nuclei is A, and the atomic number is Z. For sufficiently high values of Eq, we can neglect 
electric forces and assume that, as the nuclei approach one another, an overlap between 
the spatial distributions of the initially cold medium will appear from a certain instant of 
time onward. To obtain an approximate measure of the strong interaction which results in 
this situation, let us estimate the mean free path. The cross section for the interaction 
between the initial elementary particles, i. e., nucleons, is known from experiment (see, for 
example, [3]) and is of the order of the pion Compton wavelength h/m„c, whereas the role 
of the various constraints imposed by the Pauli principle decreases with increasing Eq. Very 
approximately, therefore, we have 

' ' (1) 



where n is the density of nucleons in nuclear matter (the possible creation of new particles 
will require a more careful analysis, and in the discussion below we shall interpret n as 
the spatial density of baryon charge). The mechanism and the possibility of a theoretical 
description of the phenomenon are very dependent on the relative free path given by (1) and 
the nuclear radius 

R = roA'/\ (2) 

It is well known that 

ro ~ h/niT^c ^ 1.4 ■ 10^"*^^ cm 

(this has already been used above in estimating n). We thus find that, in the case in which 
we are interested here^ 

l/annuR ~ A-^'^ < 1. (3) 

Since the mean free path is short, the initial stage is the fusion of the nuclei into a "compound 
system" . It is, however, important to emphasize the difference between this system and the 
usual compound nuclei formed, say, by nucleon capture. In cold or relatively low-temperature 
nuclei there is no appreciable pressure, and such nuclei exhibit no noticeable tendency to 
expand. On the other hand, a very high pressure is produced during the formation of the 
system in which we are interested. Thus, the most conservative estimates, which do not 
take into account the compression of the medium during fusion, show that the pressure is 
proportional to the total internal energy E — 2AEq. This results in the expansion of the 
compound system into vacuum. The condition given by (3) enables us to consider the second 
stage, i. e., expansion, in hydrodynamic terms. It is not clear whether any systematic theory 
would be capable of providing a detailed quantitative description of the "fusion stage" . The 
essential feature is that the entropy of the system increases from zero to some maximum 
value S. 

The formulation of the problem is thus quite close to the suggestion put forward at one 
time, on Fermi's initiative [8], for the description of collisions between relativistic strongly- 
interacting elementary particles. These ideas were extended further in an interesting paper by 



^In our previous papers [4-7] on the macroscopic treatment of apparently unrelated nuclear phenomena, 
we always came across the condition kfR ^ 1, where fc/ is the limiting momentum of the Fermi liquid 
quasiparticles. It is readily seen that the condition given by (3) reduces to a very similar criterion. 
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Landau [9]. Without going into the various aspects of this complex problem, we shall merely 
note some of the features of the difference between the nuclear case and the "elementary 
interaction" between two hadrons. The initiating interaction between the two initial particles 
may in itself be capable of creating "hadron matter" in macroscopic amounts, but this entire 
question is, to some extent, shrouded in doubt. On the other hand, during the fusion of 
heavy nuclei, the number of particles is known to be macroscopically large because of the 
nucleons that are present right from the beginning. This enhances the credibility of the 
above thermodynamic and hydrodynamic conditions. 

Strictly speaking, only the second stage of the process, i. e., the expansion stage, will be 
subjected to theoretical analysis. Let us begin with a few preliminary remarks on the physics 
of the phenomena involved in this process. The expansion of matter into vacuum occurs with 
near-sonic or ultrasonic velocities, so that viscous friction and thermal conductivity can 
hardly be expected to lead to an appreciable increase in entropy. If the resulting adiabatic 
motion {S — const) of the continuous medium is to be treated in a simplified mechanistic 
way, the internal energy of the liquid will, so to speak, play the role of potential energy. 
This provides a clear physical picture of why the overall character of the motion of the 
individual elements of the medium (the fluid particles) depends on the order of magnitude 
of the velocities communicated to them. 

In the case of inelastic collisions between nuclei, the physically interesting region is the 
relatively extensive nonrelativistic region £'0 <S m„c^ (m„ is the nucleon mass), which is 
even more accessible to current experimental possibilities. We shall write the nonrelativistic 
energy of a fluid particle in the form of the sum 

-Mv^ + e, 

i.e., the sum of kinetic and potential energies (M = const by definition and e is the internal 
energy). In very approximate calculations, we can initially ignore the energy of interaction 
with the ambient fluid (i. e., the work done by pressure), and suppose that the velocity v 
increases due to the reduction in s during the adiabatic expansion. The latter leads to a 
subsequent reduction in pressure, so that the assumption that the interaction between the 
fluid particles is small will become increasingly valid. The net result is that Mv'^ » £, i.e., 
the fluid particles become "freed" and execute inertial motion with v ~ const. This condition 
also determines the validity of the assumption that the true particles of the medium have 
negligible thermal velocities (due to cooling on expansion), as compared with the transla- 
tional velocity v of the fluid. Thus, the flnal velocity distribution of the particles, i.e., the 
reaction products, is predetermined while, on the other hand, the hydrodynamic conditions 
which demand that the mean free path is small in comparison with the linear dimensions of 
the entire system may still be valid. Essentially, this stage is actually reached in a relatively 
short time t > I/uq, where / represents the linear dimensions of the system before the onset 
of expansion, and Uq is the initial velocity of sound. 

The foregoing general ideas lose their validity as we enter the ultrarelativistic region 
Eq ^ m„c^. The single expression 

e 

(1-t;Vc2)V2 
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cannot be divided into "kinetic" and "potential" components in an entirely natural fashion. 
We note that the idea of a "freed" liquid particle is not altogether consistent with the general 
character of relativistic relationships. It is clear, for example, that the reduction in e should 
be compensated by a reduction in the denominator. The fluid continues to accelerate and, 
in reality, the pressure [for which at ultrarelativistic temperatures one usually employs the 
equation of state given by (30)] remains effective. The only process capable of terminating 
the reduction in e during expansion, and of stabilizing the velocity, is the formation of 
individual particles in the hadron matter, the rest masses of which begin to dominate all the 
contributions to the internal energy. Here again we return to the situation where the energy 
and angular distribution of the reaction products are predetermined and, correspondingly, 
the equation of state for the mcdlTim changes and departs from (30). In the opinion of 
Pomeranchuk [10] and Landau [9,11], this occurs at temperatures T ~ rriT^c^. 



2 Collisions of nonrelativistic nuclei 

We shall suppose below that the change in the internal state of the medium during expansion 
is described by the Poisson adiabatic curve [12]: 

pV^ = const. (4) 

If we recall that dE — —pdV and integrate, we can write the basic relationships in the 
following form, which is particularly convenient for subsequent calculations: 

7 = - — -j-y, w — — - — u , aw = [zv + l)uau. (5j 

In simple cases, the parameter 7 is the ratio of speciflc heats, but this is not essential; w is 
the enthalpy per unit mass. Moreover, for adiabatic (isentropic) flow 

s (X n (X p (X v?"^^, (6) 

where p is the density of the spatial mass distribution and s is the entropy per unit volume. 
For the so-called simple (self-similar) rarefaction wave, we have 

« + ^ = «o, (7) 

where u is the local velocity of sound (see, for example, [13]) and 

Vmax = (2Z^ + 1)^0 (8) 

is the limiting value of the velocity of free expansion of the medium which is initially at rest 
in vacuum. 

To obtain an estimate for the preliminary compression of nuclear matter, we shall suppose 
that the fusion of nuclei occurs gradually. Initially, in the region of space where the two media 
have come into contact, the liquid undergoes intensive "boiling" but, outside this region, it 
remains cold. Since, prior to collision, the product of the nucleon momentum by its velocity 
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is 2Eo, we find that the momentum transported out of the ambient space through unit area 
on the separation boundary per unit time is 

P = '^noEo, (9) 

where tiq is the usual equilibrium density of the baryon charge at zero temperature. The 
momentum transfer specified by (9) is obviously equivalent to a pressure p. After the medium 
has been brought to the boil, the pressure in the medium is approximately given by 

P = InEo, (10) 

which is the equation of state for an ideal gas^. Assuming that mechanical equilibrium is 
established more rapidly than thermal equilibrium in the neighborhood of the separation 
boundary, we can equate the expressions given by (9) and (10). This yields 

n/no a; 3 (11) 

prior to the onset of free expansion. We emphasize that the result given by (11) is insufficient 
to determine both the longitudinal and transverse size of the figure at the very beginning of 
the hydrodynamic stage. The medium may undergo some flow in the transverse plane which 
is perpendicular to the x axis during the fusion of the nuclei, and this is not impeded by 
external pressure. More accurate estimates of the radial size L > R reached in this direction 
are difficult because of the highly non-equilibrium character of the fusion stage (see also the 
Introduction). 

Let us now consider the adiabatic stage of the expansion process. The initial configuration 
can be schematically represented by a disk of thickness 21. It is natural to assume that 

I < L. (12) 

In the first approximation, therefore, the hydrodynamic flow can be looked upon as one- 
dimensional. The symmetry of the problem enables us to confine our attention to the 
region x > 0. In addition to the coordinate measured from the center of symmetry, it will 
occasionally be useful to use the variable x' = x ~ I. The edge of the distribution of matter 
moves forward with the velocity given by (8). As long as t < 1/uq, the situation is no different 
from the solution of the well-known problem on the expansion of a half-space into vacuum. 
Against the flow, we have the propagation of a simple wave up to the "weak discontinuity" 
x' = —Uot (i. e., the point at which the sonic signal reaches at this time; see, for example, 
[13]). When t > 1/uq, the weak discontinuity moves in the positive direction of the x axis, 
and the relative size of the region occupied by the simple wave decreases rapidly. The space 
on the other side of the weak discontinuity, where the appropriate value of the so-called 

^This means that we are neglecting the potential energy of the interaction between the nucleons. The 
assumption that an ideal gas is produced seems, at first sight, to be somewhat drastic. Nevertheless, there 
arc reasons to suppose that it does, in fact, lead to a reasonable description of the main features of the 
phenomenon. It is clear from the foregoing that the resulting particle-energy distribution essentially reflects 
the hydrodynamic character of the process, but is not too sensitive to the particular choice of the adiabatic 
curve. It is also important to remember that the contribution of the interaction energy rapidly decreases 
during the expansion process. 
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general integral of hydrodynamic equations^ is reached, begins to play the dominant role. In 
principle, a general analytic expression can be obtained for it for integral values of v [13,14]. 

In the case of an ideal gas of elementary particles, we have 7 = 5/3 and v = 1. The 
corresponding general solution can be written in the form 

xK«) = i{F.(« + ^)+ft(«-^)}. (13) 

where Fi and are some arbitrary functions. The "velocity potential" can be used for the 
implicit determination of the required functions w{x' ,t) and v{x' ,t) from the formulas 

( = 1^, .' = .|^-|^ (14) 

ow aw ov 

(see, for example, [13]). By satisfying the boundary conditions both at x' — —I {x — Qi), at 
which the fluid is at rest, and at the point of contact with the simple wave (7), we flnally 
obtain 

By substituting in (14) [see also (5)], we immediately return to the physically most interesting 
time t ^ Z/i^o (in which case, u <^U(), where uq is the initial velocity of sound in the originally 
resting medium): 

The velocity field given by (16) corresponds to the inertial motion of the fiuid particles, and 
the velocity distribution of the masses remains unaltered (see also the prehminary remarks 
in the Introduction). In fact, at any time 

pdx oc pdv oc (mq — /9)dv. (17) 

If we now transform to the new variable defined by w oc ^/e, dv oc del\fe and normalize the 
expression W{e)de oc pdx to the unit integral between and Smax^ we obtain the following 
expression for the energy distribution of the reaction products, i.e., nucleons, in the center- 
of-mass system (see Fig. 1): 

W{e)de = ^/i^(£^„, - e)^. (18) 

The presence of the cutoff point z = Zmax in (18) is a consequence of the hydrodynamic 
character of the expansion stage. The energy e averaged over the entire spectrum is given 

by 

£max = 5£ = ^Eq (19) 

•^In the model example corresponding to = 0, the "joining" of the general integral to the simple wave is 
readily achieved exactly and in an explicit form for any time t > I/uq. The fraction of energy and entropy 
which is asymptotically taken up by the simple wave turns out to be of the order of l/uot <^ 1. Similar 
estimates are characteristic for other values of i'. It must not, however, be supposed that the fact that the 
simple wave is negligible for large times f is a universal feature of all hydrodynamic problems involving the 
free expansion of material into vacuum. In the ultrarelativistic case, the fluid is rapidly accelerated and 
tends to the limiting (light) velocity so that, in general, a considerable fraction of the total energy and total 
entropy is concentrated in the simple wave. A specific example of this is the problem solved in the Appendix. 
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(e — Eq follows from energy conservation). 

When (12) is satisfied, the angular distribution of the nucleons is confined to the forward 
and backward directions. It cannot be calculated in a closed form, and we shall therefore 
confine our attention to an estimate of the characteristic angle 9 ~ Vy/v, where Vy is the 
transverse component of the fiuid-particle velocity. Its total acceleration dv/dt is given by 
the Euler equation, the transverse component of which is 

dvy ^ dp/dy fl^uj^ 

dt p ^ L ^ Lt^/^ ■ ^ ^ 

The solution given by (16) is used here for approximate purposes and is valid for U(,t < L, 
after which expansion enters the three-dimensional phase. Integrating up to the above limit, 
and recalling that v ~ i^o, we find that 

e ~ (//L)2/3. (21) 

Transforming to the laboratory system, in which one of the nuclei was at rest prior to 
collision, we obviously obtain Ei = 4Eq = As for the primary energy per nuclcon. For most 
particles, the observed angle t) between their momenta and the collision axis has the same 
order of magnitude, i.e., ■»? ~ ^ <^ 1. If we apply the Gahlean transformation to (18), we can 
readily show that 



Wier)der = [be- (ViT± (22) 



2] dei 



where ei is the laboratory nuclcon energy, and the upper and lower signs refer to particles 
travelling in the forward and backward directions in this frame, respectively. The distribution 
given by (22) is normalized to a unit total integral evaluated over both regions, and the 
corresponding branches of it are shown in Fig. 1. We note that, when £i < Ei9'^, the angular 
distribution of the nucleons becomes broad, filling the entire solid angle (in the laboratory 
system) . 

The foregoing discussion was, in fact, confined to the case Eq < 171^^0^. When the in- 
equality 

m^c^ < < m„c^ (23) 

is satisfied, the situation is modified somewhat because of the creation of a large number of 
relativistic pions during the fusion of the nuclei. In the region defined by (23), most of the 
mass is carried by the nucleons. On the other hand, the internal energy (less the nucleon rest 
mass, as is usually assumed in non-relativistic theory) resides mainly in the meson degrees 
of freedom, and these particles are also largely responsible for the pressure in the medium. 
By analogy with black-body radiation [12] [see also the next section and, in particular, the 
equation of state given by (30)], the pressure will be approximately specified by the equation 

P = InEo. (24) 

Equating the pressure given by (24) to the external pressure given by (9), which, during the 
fusion stage, describes the cold part of the medium for a certain interval of time, we find 
that 

n/no pa 6. (25) 
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The increase in the prehminary compression as compared with (11) suggests that the vahdity 
of (12) may improve ^. 

Black-body radiation and other similar ultrarelativistic modifications of matter corre- 
spond to 7 = 4/3 and v — h/2. In principle, for fractional values of there is no closed 
general analytic solution of the equations of one-dimensional hydrodynamics that are analo- 
gous to (13) and (14). However, for large times t, the asymptotic behavior of the form given 
by (16) can readily be generalized to fractional values of v. 



2 1 (26) 



r(2i/ + l) I o V- 



22-[r(i/ + 1)]2 1 



Uq - 



If we use this expression to determine the energy distribution of the particles in the center- 
of-mass system, we have for u — 5/2, 

W{e)de = ^^^% (Smax - efl'^^, Smax = 8£, (27) 

which is valid for any reaction products with nucleons and pions predominating. The quan- 
tity Smax is proportional to the mass of the particles with which we are concerned. For 
example, Sn,max/ ^■K,max = ^n/^-K- Becausc the conditions for the validity of the theory are 
unfavorable (see the last footnote), the equation = Eq is, in fact, satisfied only approxi- 
mately. Finally, if we estimate the transverse forces in the Euler equation by analogy with 
the derivation of (21) from (20) and (16), we get the expression 

e ~ (//L)i/=^ (28) 

for the effective angle at which the particles are emitted in the center-of-mass system. We 
shall not consider here the kinematics of the transformation to the laboratory system, since 
it is analogous to that discussed above for the case Eq ^ tti-,j-c . 



3 Collisions of ultrarelativistic nuclei 

When 

Eo > m„c^ (29) 

we must use the method of relativistic hydrodynamics [13,9,11]. We begin by writing 
down the basic thermodynamic relationships. In the spirit of the Landau idea [9,11] on the 
nature and the probable form of the equation of state for hadron matter at ultrarelativistic 
temperatures T 3> m^c^, we assume that 

^Nevertheless, the narrowness of the region in whieh (23) is vaUd is a serious defect of the theory applicable 
to it. The condition given by (23) may not be sufficient because the nucleon and pion rest masses are not, 
in reality, very different from one another. 
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In these expressions, e is the energy per unit proper volume of the hquid particle in its rest 
system, s is the entropy per unit proper volume, p is the pressure, and T the temperature. 
This yields the following constant value for the velocity of sound: 



u 



= c/VS. (31) 



The numerical value of k cannot be established by purely deductive means. Dimensional 
considerations suggest that 

k ~ he. (32) 

The volume of the compound system produced as a result of the fusion of the original nucleus 
is given by 

V ~ i?^^, (33) 

where R is the nuclear radius and the factor rUnC^ / Eq appears as a result of the Lorentz 
compression: 

— ~ J. (34) 

Recalling also the expression given by (2), we can readily show that the temperature and 
entropy of the system at the time preceding the onset of adiabatic expansion are given by 

T.-^,^{'I^\'\r^, S^aM'\[J^. (36) 



Let us now consider the hydrodynamic stage. Because of the geometry of the initial 
configuration, this stage has the character of one-dimensional flow over a certain interval of 
time. However, analysis shows that the increase in the influence of transverse forces gradually 
leads to the isotropization of the flow and its rapid transformation into the three-dimensional 
phase ^. As a result, the liquid is so rapidly accelerated that it becomes concentrated largely 
at finite distances from the surface which is expanding with the velocity of light. One way of 
describing this is to say that a "cavity," i.e., a region of sharply reduced density, is produced 
inside the spatial distribution of matter with this peculiar geometry. We shall first describe 
this isotropic part of the process and will return to the infiuence of the initial conditions 
later. 

It is well known that the equations of relativistic hydrodynamics are contained in the 
differential conservation laws 

3^ = 0. (36) 

where T*^ is the energy-momentum tensor of the medium [9, 11, 13, 15]. It will be convenient 
to use the system of units in which c = 1 and adopt (30) to simplify all the expressions to the 
case of spherical symmetry. Instead of the radial distance r, we shall use the independent 
variable 

^-t-r (37) 

^This phenomenon was considered qualitatively in the papers [9,11]. Landau called it "lateral" or "con- 
ical" expansion. He used conservation laws to predict a time dependence of the main quantities, which is 
confirmed by the rigorous formula (40); see below for further details. 
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and expand into a series in powers of the reciprocal of the relativistic 4- velocity 

1 



7 



»1, 



(38) 



retaining only the first two terms (i/^ ~ 7, as we shall soon show). In terms of the new 
variables, we then obtain the following set of equations for the ultrarelativistic flow: 



2'9S . ^1 2 , 2 , l-^S 1 

0' h 2— T s + 2 — s H — 



dt 



197 

7W 



0, 



Ids 3 



9f 



1 (97 



(39) 



s = 0. 



It is readily verified that all the requirements are satisfied by the following very simple 
solution^ 

The singularity at ^ = refiects, formally, the inability of matter to propagate with velocities 
in excess of the velocity of light and this is, of course, the basic feature of the equations 
of relativistic hydrodynamics. In point of fact, the value of the general integral of these 
equations given by (40) is valid only up to a certain ^0 > 0. The sphere = ^0 is a surface 
of weak discontinuity. Integrating the differential equation of its motion [it can be set up 
with the aid of the relativistic law of addition of velocities; relative to the fiuid, the weak 
discontinuity always propagates with the velocity of light which, in this case, is l/-\/3, see 
(31)], we can verify that ^ const when 



(41) 



To achieve a more specific physical interpretation of let us consider the conservation 
laws. The spatial energy density is the time component 



Tf 



00 



(42) 



of the energy-momentum tensor [9,13,15]. Moreover, the total entropy 5* is a constant in the 
case of adiabatic flow. Using (30) and (40), we obtain 



' 3 ' 3 2 



23/2 ^ 



(43) 



-^73' 

so 



^Froin the more formal point of view, the simplicity of this solution and the complexity of the one- 
dimensional Landau-Khalatnikov solution [9,11,16] are probably connected with the three-dimensional char- 
acter of real physical space. Both the hydrodynamic equations (39) and the equation of state (30), which 
is taken into account in their derivation (this also implicitly assumes the three-dimensional character of 
space), correspond to this nature of physical space. In this sense, complete concordance of the equations of 
hydrodynamics and thermodynamics in the "one-dimcnsional world" also results in an exceedingly simple 
solution. We shall not reproduce this solution here and merely note the following: when the one-dimensional 
analog of the thermodynamic relationships given by (30) is used, the equations of hydrodynamics turn out 
to be strictly linear and can be readily solved in general form. The solution is some explicit function of the 
initial conditions. 
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which enable us to calculate a-nd the arbitrary constant C: 

-S(!)''-5=S(!)'» 

Consequently [see (33), and we return to ordinary units] 

^ V'l' ^ Rirur^c^ / E,f'\ (45) 

When the temperature is reduced to T ~ rriT^c^, the individual particles are finally formed, 
and the relativistic acceleration mechanism ceases to operate (see also the preliminary re- 
marks at the end of the Introduction). To estimate the corresponding time t, let us return to 
(40). It is clear that the volume in which the medium is concentrated is ~ {ctY^Q. Moreover, 
7 ~ ct/^o- Wc may therefore conclude that the order of magnitude of the "proper volume" 
is (ct)^^o7 ~ (c^)^- III this volume, the above temperature corresponds to pion separations 
of the order of their Compton wavelength. Thus, 

ct ^ -^Nll\ (46) 

mTrC 

We note that this result is similar to the well-known formula given by (2). Since iV^ 3> A, 
comparison with (45) shows that the inequality given by (41) is clearly satisfied. The energy 
spectrum of the particles must be judged from the entropy distribution (see [9,11]). Its 
observed density is 57 and is determined by (40). Therefore, 

W{e)de oc s-^d^ = sj—dj oc 7^ci7. (47) 

This distribution cuts off sharply at 7 = •jmax, which corresponds to { = ^0 at time given by 
(46). It is readily seen that 

J^)-^ ,48) 

The particle-energy distribution in the center-of-mass system, normalized to unity, as- 
sumes the form 

3 3 

W{e)de = e^de, e = -Smax, £ < ^max = mc^lmax- (49) 

In this expression, m is the rest mass of the particular type of particles with which we are 

concerned. 

Let us now briefly consider the kinematics of the transformation to the laboratory system. 
Elementary relativistic transformation yields: 

El 
mnC^ 
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where Ei is the primary laboratory energy per nucleon in the bombarding nucleus. The 
particle energy and angular distribution can be found without great difficulty, but the process 
is laborious. Integrating it with respect to one of the variables, we find that 

p (50) 

CLEl, £l,max j-, ^maxi 
-C/Q 



W{ei)dei- ' 



1 



2 



for the angular and energy distributions, respectively. The quantity do — 27r'&d'& is the solid- 
angle element and, as can be seen, 9 ~ rrinC? / <^ 1, i.e., the angular distribution is confined 

to forward directions in the laboratory system. Wc emphasize these simple distributions are 
valid for the great majority of particles but, strictly speaking, not for all of them. As a 
matter of fact, the laboratory energy has the lower bound 

I m Eq 

^ ^max 

which is relatively low but still ultrarelativistic. In the "soft" part of the spectrum adjacent 
to Si^rnim the particles are emitted at relatively large angles right up to the maximum possible 

2 

rUnC 1^1? ^ i9 

"max — Imaxi ^ "max 

-C/Q 

These details of the ultrarelativistic distributions are illustrated in Fig. 2. 

One further remark must be introduced in connection with the foregoing. The general 
principles of solution of this kind of hydrodynamic problem would appear to enable us to say 
that the region of space < cannot be absolutely "empty." It should contain the simple 
wave which is in direct contact with vacuum. Since in the equation of state given by (30) 
the edge of the distribution of matter (strictly speaking, it, too, is a weak discontinuity) 
always moves with the velocity of light, the radial size Co of the simple wave will also remain 
constant in the ultrarelativistic limit which we have considered. To estimate it, therefore, 
we must return to an earlier stage in the expansion process. 

The one-dimensional Landau-Khalatnikov theory [9,11,16] is valid for ct ^ ,^o- The initial 
configuration was characterized by the longitudinal size / ~ RrrinC^ / Eq [see also (33)]. When 
ct > Z-\/3, both weak discontinuities move in the same, positive, direction. It is readily shown 
that, in the one-dimensional relativistic simple wave (self-similar, see, for example, [13]), we 
have ^ 

Cocx^^ A= ^^^j =7-4v^R^0.07, (51) 
which describes the distance Co between the weak discontinuities as a function of time. Thus, 

-I-R"^. (52) 
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if, for the purpose of very approximate calculations, we neglect the effect of the small ex- 
ponent, and take into account the short duration of the entire one-dimensional phase of the 
expansion process. Comparison with (45) then yields 



Co < Co- 



(53) 



We may, therefore, neglect the contributions of the energy and entropy of the simple wave, 
and this was taken into account in the derivation of the formulas considered below. 

4 Charge distribution of reaction products 

The fact that the individual particles (hadrons) have certain discrete quantum numbers, i.e., 
different "charges," enables us to derive a number of interesting relationships. 

The equilibrium character of the resulting electric charge distribution is clear even from 
(46) . Immediately after the formation of the individual hadrons, the free path ~ h/m.,^c is still 
small in comparison with the linear dimensions of the entire system, so that hydrodynamics 
and thermodynamics remain valid, as before, for an appreciable length of time even after 
transition to the region T <^ m^c^ in which we have a Boltzmann gas with a practically 
constant number of particles ^. Under these conditions, elastic interactions between the 
particles, including charge-transfer processes, are sufficiently effective. 

We shall base our analysis on the principle of isotopic invariance (sec, for example, [17- 
19]). When nuclei with the same number of protons and neutrons coalesce, the initial state 
is completely isotropic in isotopic space, with all the ensuing consequences. In particular, 
all pions (tt"*", 7r°, tt") are then created in equal numbers. However, in practice, sufficiently 
heavy nuclei have a neutron excess A — 2Z. Using the analogy with thermodynamics, and 
the statistics of rotating bodies [12], we can adhere to the point of view that, in equilibrium, 
a fluid particle rotates as a whole in isotopic space with angular velocity Vt. The Boltzmann 
distribution then contains the factor exp{/ifiT/T}, which includes the component r of the 
particle isospin along the rotation axis. Consequently, 



where Np is the number of protons and is the number of neutrons among the reaction 

products, and the other subscripts refer to pions of the appropriate type. 

The validity of relationships such as those given by (54) does not depend on the presence 
of other particles. Let us suppose now that antibaryons can be practically neglected, and 
baryons are represented only by protons and neutrons. Conservation of the baryon charge 
2A of the entire system then yields 



^Transition to the Boltzmann region T <C mTrC^ is accompanied by the strong suppression of pion anni- 
hilation processes because, as the density falls, the role of triple (and higher order) collisions falls rapidly to 
zero. 




(54) 




(55) 
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Let us now apply the conservation of electric charge Np + Ntj+ — N^- — 2Z. We have 




A-2Z A-Np _ 1 



(56) 



A + {'i/3)N^' A-2Z l + (4/3)(M^/^)' 



Thus, after the reaction, the neutron excess A — Np decreases in comparison with its original 
value A — 2Z, and hence the pion fraction contains more negative pions than positive pions. 

We note that the above formulas are even more valid for Eq < m^c^ when, roughly 
speaking, there is not enough energy for antinucleon creation. Even in the absence of pions, 
the relative neutron excess is not large enough to enable us to assume that Ml/T -C 1, as 
above. For the region defined by (23), we can readily show that, very approximately. 



In the ultrarelativistic limit Eq 3> mnC^, the situation can, at least in principle, become 
modified by the creation of baryon pairs. However, under these conditions, since N^^ ~ S 
(see (35) and [9,11]), the neutron excess in the nucleon fraction is negligible compared with 
the initial excess. 

5 Discussion 

Let us now briefly review the conclusions of the theory of collisions between energetic nuclei, 
which refer to the energies of the individual particles after interaction. Their mean value is, 
as a rule, of the order of the temperature Tq of the resulting compound system. However, 
the shape of the energy spectra of the reaction products does not in itself exclude the 
possibility that the mechanism may be interpretable as purely thermal and "evaporative." 
It is difficult to imagine, for example, that the restriction on the maximum energy of the 
emitted particle is due to anything other than the hydrodynamic character of the expansion 
of the compound system. There is a particularly sharp jump in the distribution function at 
^ = ^max in the ultrarelativistic limit (see (49) and the explanation in text). When we refer 
to the nonrelativistic case Eq <^ fnnC^, we must also emphasize the shape dej ^ of the soft 
part of the spectrum, which is totally uncharacteristic for particle-evaporation processes in 
the case of the usual compound nucleus. When the necessary experimental data become 
available, therefore, one would hope to be able to achieve a sufficiently reliable identification 
of the hydrodynamic mechanism discussed in the present paper. 

We must now briefly consider the spcciflc features of collisions that are not of the head-on 
type. The compound system whose evolution is described by the above theory arises in the 
region of space where the coUiding nuclei overlap. Those parts of the nuclei which do not 
overlap remain as relatively cold fragments, in effect, truncated on collision. They largely 
continue to exccTitc incrtial motion with energy Eq per nucleon. Subsequently, the shape 
of a fragment in its rest system tends to an equilibrium, and the oscillations of the surface 
become transformed into heat. The flnal temperature reached in the course of this process 
is probably a slowly- varying function of the primary energy Eq and is low. Consequently, 
the velocity of the nucleons evaporated from the fragment is also small in comparison with 




(57) 
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its translational velocity as a whole. The nucleons evaporated by this mechanism should 
therefore produce an additional monochromatic peak at e ~ i?o in the energy spectrum (in 
the center-of-mass system). This interesting feature of the phenomenon suggests that the 
experimental energy distributions should be even more informative. 

We would hke to express our gratitude to A.I. Baz', I.I. Gurevich, L.P. Kudrin, V.A. 
Novikov, A.A. Ogloblin, I.I. Roizen, Ya.A. Smorodinskii, Yu.A. Tarasov, and K.A. Ter- 
Martirosyan for discussions of the present results. 



Appendix 

There is some methodical interest in the problem of expansion of matter which initially occu- 
pies uniformly a spherical volume of radius R at rest. We shall assume that the temperature 
is ultrarelativistic. 

In the limit t ^ R, when the one-sided expansion away from an internal weak discon- 
tinuity = .^0 has taken place, we have the general integral given by (40). On the other 
side, C < Co we have a spherically symmetric simple wave. To establish the shape of the 
singularity on the surface of the external weak discontinuity, i. e., on the boundary with 
vacuum, let us consider the corresponding self-similar solution (which depends only on the 
variable r] = r/t). Detailed analysis, which we shall omit for lack of space, leads to the 
following natural-looking result 



Here, in contrast to ^, the coordinate ( is measured from a different bounding surface of the 
light signal (see also below); Co is the position of the internal weak discontinuity on the C- 
scale. In (A.l). the constants have been chosen so as to ensure that it agrees with (40) when 
^ = ^0- Foi' the conserved total energy and entropy, we now have E = E' + E", S = S' + S", 
where E' and S' are given by (43) and E" and S" can be calculated by analogy, using the 
solution (A.l) for the simple wave, and then integrating with respect to ( between zero and 
Co- The result is 



To ensure that the "asymptotic" solution (40), in general, misses the region in which the 
equations of hydrodynamics are at least formally satisfied (this is discussed below with a 
suitable choice of the origin of time), the sphere ^ = must be reduced to a point. It is clear, 
on the other hand, that, in reality, this kind of singularity at r = occurs only when the 
surface of the internal weak discontinuity contracts to the origin, t — Ry/S. During the same 
time, the external weak discontinuity will move forward to a distance RV^. Subsequently, 
the distance 

C-e = Co-eo = (V3-l)i? (A.3) 

between the singular surfaces (the spheres C = and ,^ = 0) of the two solutions to be 
matched undergoes no changes because both surfaces expand with velocities strictly equal 
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to the velocity of light. This is the necessary third condition which, together with (A. 2), 
determines finally all three arbitrary constants C, ^q, and (q. Since 

we obtain an algebraic equation of a high (sixth) degree in the ratio R/^o- Numerical solution 
yields 

Co = 0.52/?, Co = 3.25i?, Co/Co = 6.3, 

24x/2(l + 5Co/3eo)^ ^3c, 
1257r (1 + Co/eo)^ 

Returning now to ordinary units, and using the dimensionless small combinations 

, ^ r + V^R „ ^ r-R 

z = 1 , z = 1 

ct ct 

for the sake of brevity, we finally obtain 

t > R/c. 

Equations (A. 4) and (A. 5) form the solution of our problem. 

If, as a result of further expansion, the medium cools down to nonrelativistic tempera- 
tures, the equation of state given by (30) will be violated and the particle energy distribution 
will cease to vary. We shall not consider the details of this; the necessary derivations are 
similar to those leading to (47) and (49). Wc merely note that the "weakness" of the internal 
discontinuity at ,^ = is reflected in the continuity of the thermodynamic and hydrodynamic 
quantities (but not of their spatial derivatives). However, to determine the particle-energy 
spectrum, we must transform to the 7 scale, in which case the spectral density W{e) itself 
exhibits a discontinuity. It is not difficult to show that 

directly at the singularity e = Sq = mc^7(,^o) (the subscripts + and — indicate the values of 
the functions to the right and to the left of it). The final result is that the energy density 
W{e) increases discontinuously [see (A. 4)] and thereafter decreases in accordance with the 
formula W{e)de oc de/e^, e > Sq. On the other hand, it follows from (53) that, in the case of 
the problem considered in Sec. 3, the initial geometry ensures that the distribution function 
W{e) falls to a negligible value for e = eQ = Smax- 
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Figures captions 



Fig. 1. The dashed curve shows the particle spectrum in the c.m.s. Curve 1 refers to 
particles travelling in forward directions in the laboratory system. Curve 2 refers to particles 
travelling in the backward directions. In the laboratory system eimax = (6 ± 2y/5)e. 

Fig. 2. These graphs were plotted for (48) replaced by an equation in which the propor- 
tionality factor was taken to be equal to unity. 
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This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/0312102vl 



This figure "fig2.gif" is available in "gif" format from: 
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